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O (57) Abstract: The present invention provides a predictive tool for predicting Apnea Hypopnea Index (AHI) in the diagnosis of 
^ Obstructive sleep apnea (OS A). The predictive tool is developed by recording pulse oximetry readings (104), and obtaining delta 
Q index, oxygen saturation times and oxygen desaturation events from the oximetry readings (106). A multivariate non-parametric 
analysis (110) and bootstrap aggregation is performed to obtain predictive models (112) which can be used to predict AHI in an 
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METHOD FOR PREDICTING APNE A-HYPOPNEA INDEX FROM 
OVERNIGHT PULSE OXIMETRY READINGS 

5 

This application claims priority to U.S. Provisional Application No. 
60/505,288, filed September 23, 2003, the disclosure of which is incorporated herein 
by reference. 

10 FIELD OF THE INVENTION 

The present invention relates generally to the field of sleeping disordered 
breathing. More particularly, the present invention provides a predicted tool to 
predict the apnea-hypopnea index. 

15 DISCUSSION OF RELATED ART 

The obstructive sleep apnea (OS A) syndrome is a major health problem 
affecting 2-4% of the middle-aged population. At present, polysomnography (PSG) 
is considered the reference standard diagnostic test for this condition (AHCPR 
Publication 99-E002, February 1999). However, PSG is costly and time consuming. 

20 As a result, primary care providers may be reluctant about ordering PSG and 

patients unwilling to attend their tests. Overnight pulse oximetry has been proposed 
as a simpler alternative to PSG in the diagnosis of OSA because it is readily 
available, relatively inexpensive, and could potentially meet the large demand for 
diagnostic testing in the community. It can be easily done at home and repeated, if 

25 need be, which is not the case with a PSG even performed at home. 

At present there is no definite established AHI criterion for the diagnosis of 
OSA with the threshold varying from 5-20 events/hour. Most clinicians will modify 
initiation of treatment depending on the patient's symptoms and other clinical 
characteristics. Recent reports suggest that even what is considered as mild sleep- 

30 disordered breathing is associated with hypertension and cardiovascular disease, and 
that these patients with mild disease may also benefit from continuous positive 
airway pressure (CPAP) therapy. A consensus statement recommended that 
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treatment be administered if the AHI is >30 events/hour regardless of symptoms. 
However, results from a recently published study do not support this 
recommendation since patients with an AHI of >30/hour who did not have daytime 
sleepiness did not benefit from CPAP therapy. Therefore, a prediction of the actual 
5 AHI from overnight oximetry would be more clinically meaningful than a 

dichotomous answer to the presence of OSA, and that it would be more useful if this 
prediction can be computerized to eliminate the problem of interobserver and 
intraobserver variability. 

Several quantitative indices derived from overnight pulse oximetry have 

10 been used to predict the presence of OSA. These indices include the number of 
oxyhemoglobin desaturation events below a certain threshold, usually 3% or 4% 
decline from baseline, the cumulative time spent below an oxyhemoglobin 
saturation of 90% (CT 90), and the delta index, a measure of the variability of the 
oxyhemoglobin saturation (Levy et al, Chest (1996); 109:395-99). One study has 

15 suggested that the number of desaturation events >4% as well as the 12-second delta 
index also predicts the response to CPAP therapy in patients with OSA (Choi et al., 

r 

Respir Med (2000); 94:895-99). Although these quantitative indices appear to hold 
more promise than visual inspection of the overnight pulse oximetry tracing, there 
has been no systematic comparison of their relative utility in the diagnosis of OSA. 
20 As a result, physicians select different parameters to interpret overnight pulse 
oximetry results. 

Most published studies utilizing these quantitative oximetry indices have 
been performed at a single institution. Thus, the applicability of these indices to the 
general population remains uncertain. In addition, their accuracy has been validated 

25 using different threshold values of the AHI due to a lack of established criterion for 
the diagnosis of OSA. In practice, most physicians tend to modify the initiation of 
treatment for OSA depending on the patient's symptoms and clinical characteristics. 
Prediction of the actual AHI from overnight pulse oximetry is likely to be more 
useful than using threshold values to define OSA that has been customary in all but 

30 a few studies involving pulse oximetry. 
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SUMMARY OF THE INVENTION 

The present method provides a method for predicting Apnea Hypopnea Index 
(AHI) using the predictive tool described herein. Based on the predicted AHI, an 
individual can be diagnosed as having or not having obstructive sleep apnea. 
5 The predictive tool is obtained by identifying patients having obstructive 

sleep apnea; determining actual AHI in the patients using standard methods; 
obtaining pulse oximetry recordings from the patients; determining a set of indices 
from the pulse oximetry readings such as delta index, time spent at selected oxygen 
saturation levels, and oxygen desaturation events; performing random sampling with 

10 replacement to generate different data sets; performing multivariate non-parametric 
analysis in the data sets to generate prediction models; obtaining a predictive AHI 
value for each patient corresponding to each prediction model; and comparing the 
predictive AHI value to actual AHI value for each patient to provide weight to each 
prediction model. The prediction models and the weight assigned to each model 

15 provides a predictive tool for determining AHI in an individual. 

In another embodiment, a method is provided for determining AHI in an 
individual. For the determiniation of the AHI, pulse oximetry readings are obtained 
from the individual, and a set of indices determined comprising delta index, time 
spent at selected oxygen saturation levels and oxygen desaturation events. This data 

20 is then input into the prediction models and based on the weight assigned to each 

prediction model a weighted average (final value) is determined for a AHI value for 
that individual. Based on this final value of AHI, a classification of the individual 
as having or not having OS A can be made. 

In another embodiment, a software storage device is provided having thereon 

25 a computer readable code for executing steps to determine the AHI in an individual 
as described herein. 

In another embodiment, the predictive tool further utilizes a set of indices for 
generating predictive models selected from the group consisting of data obtained by 
standard overnight PSG with recordings of electroencephalogram (EEG), electro- 

30 oculogram, submental and bilateral leg electromyograms, electrocardiogram, airflow 
measured qualitatively by an oral-nasal thermistor and respiratory effort measured 
by thoracoabdominal piezoelectric belts. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a representation of the steps for developing the predictive tool of 
the present invention. 
5 Figure 2 is a representation of steps involved in classifying a patient as 

having OSA or not having OSA. 

Figure 3 is a representation of the aggregated model using a combination of 
all ten oximetry indices to predict the AHI from the derivation group (n=224). The 
predicted and actual AHI values are shown on logarithmic scale. The coefficient of 
10 determination (r 2 ) between the actual and predicted AHI was 0.70 which is 

significantly improved compared to 0.60 (data not shown) using the initial model 
(p<0.05). 

Figure 4 is a representation of the positive likelihood ratios (+ SE) in: A) 
Derivation Group, B) Validation Group I, and C) Validation Group II according to 
15 the severity of the AHI derived from multivariate prediction models using the delta 
index alone (filled triangles) versus the aggregated model (filled circles). The 
aggregated model improved the prediction mainly at both ends of the AHI spectrum 
(<5/hour and >30/hour). 



20 DETAILED DESCRIPTION OF THE INVENTION 

The following abbreviations have been used in this application. 



25 



30 



AHI 


apnea-hypopnea index 


ASC 


Associated Sleep Center 


BMI 


body mass index 


CI 


confidence intervals 


CPAP 


continuous positive airway pressure 


CT 


cumulative time 


EEG 


electroencephalogram 


MARS 


multivariate adaptive regression splines 


OSA 


obstructive sleep apnea 


ROC 


receiver operator characteristic 


SE 


standard error 
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VAMC Veterans Affairs Medical Center 

References relevant to the invention are listed at the end of the specification. 

The term Apopnea Hypopnea Index or AHI as used herein means as the 
5 number of apneas plus the number of hypopneas per hour. 

The term delta index as used herein means the absolute deviation from the 
mean of the oxygen saturation measured at twelve second intervals. 

The present invention provides a predictive tool for determining AHI from 
overnight pulse oximetry readings. From the AHI it can be diagnosed whether or 
10 not an individual has OS A. The development of the tool is based on a study in 

which it was observed that: 1) among the different oximetry indices, the delta index 
was the best predictor of the presence of OSA, although desaturation events 
provided similar levels of diagnostic accuracy, 2) the delta index had good 
sensitivity but low specificity, 3) a bootstrap aggregation of models involving a 
15 combination of all the oximetry indices (compared to using the delta index alone) 
improved the precision of the prediction of the AHI, and 4) the prediction model 
developed in this invention was validated in two independent sleep clinics. To our 
knowledge, there has been no previous study that has compared systematically the 
relative utility of the various quantitative indices derived from overnight oximetry in 
20 the diagnosis of OSA. 

To prepare the predictive tool of the present invention (Figure 1), patients 
with OSA are identified (Step 100). For each individual, polysomnography (PSG) 
recordings are carried out, preferably overnight. Standard overnight PSG recordings 
include electroencephalogram (EEG) recordings, electro-oculogram and 
25 electromyogram (EMG). Airflow recordings and respiratory effort can also be 

recorded. During these recordings, sleep stages are scored in selected time epoch. 
Each epoch is analyzed and AHI is determined by standard parameters (Step 102). 
In an alternative embodiment, electrocardiogram readings may also be recorded. 

During the PSG recordings, pulse oximetry readings are also obtained. For 
30 these recordings, a patient is connected to a pulse oximeter and overnight readings 
are obtained (Step 104). The pulse oximeter may be a portable unit. The oximeter 
probe is applied to the patient's finger overnight. The device can be turned off in the 
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morning and the data mailed back to the health care facility for analysis. The data is 
downloaded from the pulse oximeter into a computer. Typically this is done in an 
ASCII format. Then the data is analyzed as described herein. The data is first 
checked for artifacts as described by Taha and colleagues (Taha et al., Sleep (1997), 
5 20:99 1 - 1 00 1 ) and then the total recording time, the delta index, the cumulative time 
spent under an oxygen saturation of selected percentages and the number of 
desaturation events at selected percentages are determined (Step 106). In one 
embodiment, the selected percentages for determining time spent under oxygen 
saturation are generally below a certain levels of oxygen saturation between the 

10 range of 50 to 100%. For example, time spent under oxygen saturation of 80, 81, 
82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99 and 100% can be 
determined. In one embodiment, the time spent under oxygen saturation of 80, 82, 
84, 86, 88 and 90% is determined. For desaturation events, the level of desaturation 
may be selected from 2 to 20%. For example, in one embodiment, the level of 

15 desaturation is 2, 3 and 4%. 

A sufficient number of random data sets are generated from the original data. 
Each of the random data sets utilize the process known as bootstrapping. This 
process allow for selection of random patients to fill a set of equal size to the 
original data set. However, each time a selection of a patient is made, the patient is 

20 then put back into the pool of available patients for the selection of the next patient. 
Therefore, it is possible the same patient may appear in the same data set multiple 
times or none at all. The ability to be able to derive many random sets from a 
limited number of available patients affords the mathematical flexibility necessary in 
the development of the model (Step 108). In one embodiment, the number of data 

25 sets generated in 20. 

The data are then subjected to multivariate non-parametic analysis to 
determine relationships of indices in each data set (Step 110). This can be done by 
standard statistical methods or by using commercially available software (such as 
Multivariate Adaptive Regression Splines, Salford Systems, San Diego, CA). 

30 Multivariate adaptive regression splines (MARS) were used to develop prediction 
models. The splines consist of one or more of a series of linear segments joined at 
adjacent ends that could be fitted to nonlinear data. MARS is a multivariate 
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nonparametric procedure that builds flexible regression like models using exhaustive 
search techniques to test necessity of different predictors. Interactions between 
independent variables are simultaneously tested. From these relationships, a 
prediction model is developed for each set (Step 112). The model is adaptive 
5 because it overfits the data, and then determines the size of the model that optimizes 
the tradeoff between accuracy (bias) and variance (precision) using a ten-fold cross- 
validation. The final model is obtained through backward elimination to the optimal 
model size. Predicted value is derived as linear combination of basic functions. 

The various indices of overnight oximetry are correlated with each other. As 

10 a result, a model developed on the basis of the interrelationship of the indices is 
likely to depend upon the individual cases used to develop the model. To address 
this issue, a technique that uses an aggregated result from multiple models can be 
used (Breiman, Ann Stat (1996); 24:2350-83). A provisional predicted AHI value is 
obtained for each patient corresponding to each model and the provisional predicted 

15 AHI value is compared to the actual AHI value (determined by standard methods 
from the PSG recordings) to provide weight for the particular data set model (Step 
114). The final predictive values can lie between 0 and 360/hr beyond which they 
are truncated at the midpoint value since values outside this range are unachievable. 
Based on the weight for each model, a weighted average for predictive values (the 

20 final AHI value) is obtained for each patient (Step 116). The plurality of predictive 
models as described herein (bootstrapping) and obtaining a weighted average of the 
predicted AHI values (bagging) provides a predictive tool for predicting AHI values 
in individuals. The prediction model can be validated by obtaining data from 
another group of patients or patients from another facility. 

25 The aggregation method used herein results in a significant improvement of 

the precision of the predicted AHI. This improvement is important because even 
small improvements in precision can increase the confidence in the prediction. An 
advantage of the present method includes that it is potentially useful to identify 
severe OS A or confirm the absence of OS A; the instrumentation required for it is 

30 not expensive; and it can be performed in a patient's home. 

Other advantages are that oximetry is a simple procedure so that a repeat test 
in the patient's home on a separate night in the event of technical difficulties would 
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be reasonable. Portable oximeters capable of storing data over a prolonged 
recording are now readily available. The oximeter can be sent home with the patient 
and after the overnight recording, mailed back to the sleep laboratory or physician's 
office for downloading of the data and a computer generated report of the predicted 
5 AHI with its 95% confidence interval reported back to the physician through an 
automated system. 

The clinical utility of pulse oximetry can be assessed quantitatively from the 
likelihood ratios. The Bayesian approach is to multiply the pre-test odds by the 
likelihood ratio to determine the post-test odds. Positive likelihood ratios that cause 

10 large changes in the likelihood of disease are >10, moderate changes 5 to <10, small 
changes 2 to <5 and trivial changes <2. The data presented in Example 1 show that 
pulse oximetry using the aggregated method can produce larger shifts in the pre-test 
to post-test certainty of predicting the probability of OS A in the normal and severe 
ranges of disease severity as compared to the mild to moderate range. 

1 5 The present invention is also directed to a computational system including 

software storage devices such as a floppy disk, hard drive, a CD, a zip drive, flash 
card, smart drive etc. having thereon computer readable code for causing a computer 
to execute all or a substantial portion of diagnostic method using a predictive tool as 
described herein. 

20 The present invention also provides a method of using the predictive tool for 

determination of AHI in a patient (Figure 2). The method comprises obtaining pulse 
oximetry data from an individual (Step 200), preferably overnight, and calculating 
the delta index, time of O2 saturation at selected percentages, and O2 desaturation 
events (Step 202). In one embodiment, additional indices selected from the group 

25 consisting of: electroencephalogram (EEG), electro-oculogram, submental and 

bilateral leg electromyograms, electrocardiogram, airflow measured by an oral-nasal 
thermistor, respiratory effort measured by thoracoabdominal piezoelectric belts and 
responses to sleep history questionnaires may are obtained. These indices are used 
as input parameters for each prediction model of Step 1 12 in Figure 1 to obtain a 

30 provisional AHI index (Step 204) for each model. A final AHI value is determined 
by averaging the weighted values from each model (Step 206). Based on the final 
predicted AHI value, an individual can be classified as having OSA or not having 
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OSA (Step 208). 

The invention is further described by the examples provided herein which 
are intended to be illustrative and should not be construed as being restrictive in any 
way. 

5 

Example 1 

Five hundred and sixteen patients suspected of having OSA were enrolled 
into this study. Patients were recruited from two independent sleep clinics in 
Buffalo, New York: the Associated Sleep Center (ASC) and the Buffalo Veterans 

10 Affairs Medical Center (VAMC) Sleep Center. The eligibility criteria were all 

patients who underwent an overnight PSG for suspected sleep apnea The exclusion 
criteria were age <18 years, the sleep study was done while on oxygen 
supplementation, or CPAP titration was performed on the same night as the 
diagnostic study (split-night study). Of the 5 16 patients, 224 were assigned to the 

15 Derivation Group, 101 to Validation Group I and 191 to Validation Group n. The 
patient characteristics of the Derivation and two Validation Groups are shown in 
Table 1. 

Table 1 

20 Summary of patient characteristics. 



Patient characteristics 


Derivation 
Group 
n=224 


Validation 
Group I 
n=101 


Validation 
Group II 
n= 191 


Age (years) 


48.9 ± 12.3 


51.8 ±11.5 


56.0±12.8 §§ 


BMI (kg/m z ) 


32.3 + 7.4 


32.7 ±7.1 


32.9 ± 8.7 


Neck circumference (cm) 


41.6 ±4.5 


42.6 ± 4.8 


43.7 ± 4.2 s 


AHI (events/hour) 


18.2 ±20.0 


20.2 ±19.5 


18.2±21.2 


AHI (events/hour) range 


0-108.6 


0-92.1 


0-147.2 


AHI <5/hour 


33.0% 


19.8% 


25.7% 


AHI 5 to <15/hour 


23.2% 


33.7% 


35.6% 


AHI15to<30/hour 


25.5% 


20.8% 


18.3% 


AHI >30/hour 


18.3% 


25.7% 


20.4% 
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Definition of abbreviations: BMI = body mass index; AHI = apnea-hypopnea index 
Values represent mean + standard deviation 

■ 

§§ significantly different from derivation group and validation group I (p<0.05) 
§ significantly different from derivation group (p<0.05) 

5 

All groups had similar body mass index (BMI) and AHI. The patients in 
Validation Group II were significantly older compared to the Derivation Group and 
Validation Group I, and had a larger neck circumference compared to the Derivation 
Group. 

1 0 All patients underwent standard overnight PSG with recordings of 

electroencephalogram (EEG), electro-oculogram, submental and bilateral leg 
electromyograms, and electrocardiogram. Airflow was measured qualitatively by an 
oral-nasal thermistor and respiratory effort by thoracoabdominal piezoelectric belts. 
All signals were collected and digitized on a computerized PSG system (ASC: 

15 Rembrandt, Aerosep Corp., Buffalo, NY; Buffalo VAMC: Acquitron, Mallinckrodt, 
St. Louis, MO). 

Sleep stages were scored in 30-sec epochs using standard criteria. Each 
epoch was analyzed for the number of apneas, hypopneas, EEG arousals, 
oxyhemoglobin desaturation, and disturbances in cardiac rate and rhythm. Apnea 

20 was defined as the absence of airflow for at least 1 0 seconds. Hypopnea was 

defined as a visible reduction in airflow lasting at least 10 seconds associated with 
either a 4% decrease in arterial oxyhemoglobin saturation or an EEG arousal. An 
arousal was defined according to the criteria proposed by the Atlas Task Force. 
Apneas and hypopneas were classified as obstructive if respiratory effort was 

25 present and central, if respiratory effort was absent during the event. The AHI was 
defined as the number of apneas and hypopneas per hour of sleep. Only one person 
in each sleep laboratory, blinded to the offline analysis of pulse oximetry data, 
scored the sleep studies. 

Pulse oximetry data was collected as part of the PSG. Measurement of 

30 arterial oxyhemoglobin saturation was performed with a pulse oximeter (ASC: 
Nellcor N-200, Nellcor Puritan Bennett, St. Louis, MO; Buffalo VAMC: Biox 
3740, Ohmeda, Boulder, CO). The oximeter probe was placed on the patient's 

-10- 
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finger. Oximetry data was digitized and collected at 8 and 10Hz at the ASC and 
Buffalo VAMC, respectively, into a computerized PSG system along with the other 
sleep study parameters. The recording time was defined as lights-off to lights-on 
(approximately 10 p.m. to 6 a.m.). Recording time was used as the denominator for 
5 the various indices of overnight pulse oximetry rather than total sleep time since 
EEG is not be available when oximetry is to be utilized outside of the sleep 
laboratory. The oximeters in both sleep laboratories employed a moving average of 
3 seconds. The oximetry data was then extracted from the computerized PSG 
system for further off line analysis. Oximetry data was averaged over two-second 

10 sampling intervals. Artifacts were removed by eliminating all changes of oxygen 
saturation between consecutive sampling intervals of greater than 4% per second, 
and any oxygen saturation less than 20%. The definition of a desaturation event was 
also based on the work of Taha and colleagues (Taha et al., Automated detection and 
classification of sleep disordered breathing from conventional polysomnography. 

15 Sleep 1997 20:991-1001) . Every data point was examined sequentially to 

determine if criteria were met to define an event. The criteria for an event were a 
decrease of at least the set amount (2%, 3%, or 4%) in oxygen saturation from the 
initial data value for at least 10 sec, and at a rate greater than 0.1% sec. In addition, 
the oxygen saturation returned within 60 sec to within 1% of the initial value, or 

20 increase from its nadir by 1 .5 times or more of the set amount of the dip. Once the 
criteria were met, a new search for an event was initiated at the next data point after 
the event. The delta index was calculated as the average of absolute differences of 
oxygen saturation between successive 12-sec intervals (sum of the absolute 
differences between two successive points, divided by the number of intervals 

25 measured). The index quantifies oxygen saturation variability. In OSA, oxygen 
desaturation events associated with respiratory events cause fluctuations in the 
oxygen saturation signal leading to high delta index values. A total of 10 indices 
were calculated for each patient in the Derivation Group. The delta index, number 
of desaturation events (to 2%, 3% and 4% levels) per hour of recording time, and the 

30 cumulative time (CT) spent below 90, 88, 86, 84, 82, and 80% saturation as a 
proportion of total recording time. 

The various indices of overnight pulse oximetry (predictor variables) were 

- 11 - 
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subjected to multivariate non-parametric analysis. The prediction equation using the 
delta index alone derived using the multivariate model in the Derivation Group is: 

log 10 (AHI + 1) = 1.306 + 0.269*BF1 - 2.316*BF2 (equation 1) 

5 

where BF1 is the maximal value of either zero or (DX - 0.57), BF2 is the maximal 
value of either zero or (0.570 - DX), and DX is the delta index. The coefficient of 
determination (r^) between the actual and predicted AHI was 0.60 using this initial 

1 0 multivariate prediction model involving the delta index alone. 

To minimize the effect of differences in the correlation of various indices in 
different individuals, the predictions of different plausible models (representing 
different interrelationships between indices) were averaged as described above. 
Bootstrap aggregation model averaging ("Bagging") was used by developing 20 

15 random samples (with replacement) from the original data set. Each of the twenty 
data sets has the same size as the original derivation data set. Because random 
sampling with replacement was used, a particular patient could occur more than 
once in any of the twenty generated data sets, and some may not appear at all. For 
each of the 20 data sets, a MARS model was generated in a similar form to that 

20 shown in equation 1 . 

For every patient in original derivation data set (subscript i), and each of 
predictive models (subscript m), the predicted value of AHI 0 < pi >m < 360, was 
determined, so that every patient was assigned 20 predictions. Predictions beyond 
this range were truncated at the endpoint values because results outside this range 

25 are unachievable. The maximal value 360/hr would indicate continuous apnea since 
apnea is defined as an event of 10 sec or more in duration. 

The multiple linear regression model with bootstrap sample predictions as 
independent variables was fitted to transformed response: 

30 logio(AHI + 1) = bo + bipi . . . + b2oP20 + e (equation 2) 

wherein AHI was measured for that patient by overnight polysomnography and e is 
the error term. The best regression model was found with all subsets method. If the 
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model was not included in the final regression model, the corresponding coefficient 
was assigned to be zero. Weighted average of the predictions was used to obtain a 
single aggregated prediction of AHI for a particular patient: 

5 logio(AHI + 1) = E(b m * p m ) / E b m (equation 3) 

The aggregated model using a combination of the various oximetry indices 
increased the r* between the actual and predicted AHI (r 2 = 0.70) which was 
significantly higher than the initial model (p<0.05) as derived by equation 1 . The 

1 0 predicted and actual AHI of the 224 patients in the Derivation Group using the 

aggregated model are shown in Figure 3. In comparison with the derivation model 
using only the delta index, there was an improvement in terms of diagnostic 
accuracy using >15/hr to define OS A with the aggregated model. The area under 
the ROC curve was increased to 0.9 ± 0.02 with a sensitivity of 90% (CI: 82-95%) 

1 5 and a specificity of 70% (CI: 62-78%) using the aggregated model. The term 

"sensitivity" as used herein refers to the true positive rate and the term "specificity" 
as used herein refers to the true negative rate. 

The prediction model developed from the Derivation Group as described 
herein was validated in two independent facilities (Validation Group I and 

20 Validation Group II). The characteristics of the patients in the Derivation and two 
Validation Groups were compared using nonparametric one-way analysis of 
variance. If a significant difference was found, a multiple comparison procedure 
(Dunn's method) was used to determine a source of the difference (Sigmastat, SPSS 
Inc., San Rafael, CA). In Validation Group I, the actual AHI values of 92 out of 101 

25 patients were within the CI of the predictions for the AHI using the aggregated 

model. The proportion of patients within the confidence interval of the prediction 
was 90% (CI: 83-95 %). In Validation Group H, actual AHI values of 174 out of 
191 patients were within the CI of the predicted AHI. The proportion of patients 
within the confidence interval of the prediction was 91% (CI: 86-95 %) . 

30 To determine exactly at what levels of disease severity the delta index model 

and the aggregated model differed, we stratified the data into four groups according 
to the AHI measured by PSG: normal (AHI<5/hr), mild (5-<15/hr), moderate (15- 
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<30/hr) and severe (>30/hr). Figure 4 shows the positive likelihood ratios in all 
patient groups according to the severity of the AHI of the multivariate prediction 
models using the delta index alone versus the aggregated model. Both the delta 
index and the aggregated models improved the prediction mainly at both ends of the 
5 AHI spectrum (<5/hour and > 30/hour). The aggregated model was superior to the 
delta index model in the severe level in the Derivation Group and the severe level 
and normal level in Validation Group I, but no different in Validation Group II 
(Figure 4). 

The likelihood ratios are computed from the post test probability divided by 

10 the pre test probability and produce an indicator to assist the physician in 

determining the presence or absence of OSA in an individual patient. The likelihood 
ratios for the aggregated model in the Derivation Group were 6.9 (95%CI: 4.4 - 
10.8), 4.3 (95%CI: 2.8-6.5), 3.6 (95%CI: S.5-5.4) and 12 (95%CI: 7.1-20) for 
normal, mild, moderate and severe disease severity respectively. The likelihood 

1 5 ratios fell within the 95%CI for the normal and mild levels in Validation Group I 
and for the moderate and severe levels in Validation Group II (Figure 4). The 
likelihood ratios for the delta index model in the Derivation Group were 5.4 
(95%CI: 2.2-12.6), 2.5 (95%CI: 1.7-5.3), 3.0 (95%CI: 2.1-4.4) and 4.0 (95%CI: 
3.6-6:5) for normal, mild, moderate and severe disease severity respectively. The 

20 likelihood ratios fell within the 95%CI for the mild level only in Validation Group I 
and for all levels in Validation Group II (Figure 4). 

To obtain an overall estimate, we combined the results from the Derivation 
and both Validation Groups for the aggregated model. The overall likelihood ratios 
for the aggregated model in the Derivation Group were 4.2 (95%CI: 3.3 -15.3), 3.4 

25 (95°/oCI: 2.7-4.3), 3.0 (95%CI: 2.2-4.1) and 6.7 (95%CI: 4.9-9.2) for normal, mild, 
moderate and severe disease severity respectively. The overall likelihood ratios for 
the delta index were 3.3, 2.5, 3.0 and 4.9/hr for normal, mild, moderate and severe 
disease severity respectively. Therefore, the likelihood ratios were at the lower 95% 
CI in both the normal and severe levels of disease severity. 

30 A recent study using an automated analysis of oximetry data and a 

desaturation event definition of >4% lower than baseline, reported a very high 
sensitivity of 98% and specificity of 88%. However, this study used a definition of 
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arousals that differs substantially from the criteria proposed by the Atlas Task Force 
that has come into general use in the United States. As a result, their definition of 
hypopnea will differ substantially from ours. These investigators found that the 
addition of arousal-based scoring criteria (using their definition of arousal) for 
5 hypopnea causes only small changes in the AHI. However, a large study has found 
that incorporating arousals based on the Atlas Task Force criteria on the hypopnea 
definition do impact on the value of the AHI. Table 2 shows the comparison of our 
results to others using an AHI cut-off value of >15 events/hour to define the 
presence of OS A. Our results are consistent with others in the field although our 

10 specificity was higher using the aggregated model compared to the previously 
published studies using the delta index. 

The study by Levy and colleagues reported that the correlation between the 
delta index and actual AHI was 0.72, whereas in Olson's study the Spearman's 
correlation coefficient between the delta index and actual AHI was 0.71 . In our 

1 5 study, the correlation (expressed as Pearson's correlation) between the predicted and 
actual AHI was 0.77, which improved to 0.83 when we used a combination of the 
oximetry indices. Therefore, our prediction model provides an improvement 
compared to using a simple regression between the delta index alone and actual 
AHI. 

20 

Table 2 

Comparison to the results of others. 



Study 


Method used 


Sensitivity % 
(CI)* 


Specificity % 
(CD* 


Vazquez et. al. 


4% desaturation 


98 


88 


Levy et. al. 


delta index >0.6 


98 (96-100) 


46 (37-55) 


Olson et. al. 


delta index >0.4 


88 


40 


Present study 


4% desaturation 


94 (87-98) 


44 (35-52) 


delta index >0.63 


91 (84-95) 


59 (49-69) 


AHI predicted 
from aggregated 
model 


90 (82-95) 


70 (62-78) 



'based on an AHI cut off value >15 events/hour 

-15- 



WO 2005/030048 

CI = 95% confidence interval 



PCT/US2004/031178 



While this invention has been described through specific embodiments, 
routine modifications will be apparent to those skilled in the art and such 
5 modifications are intended to be within the scope of the present invention. 
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Claims: 

1 . A method for developing a predictive tool for predicting Apopnea Hypopnea 
Index (AHI) in an individual comprising the steps of: 

a) identifying patients having obstructive sleep apnea; 
5 b) determining actual AHI in the patients; 

c) obtaining pulse oximetry recordings from the patients; 

d) determining patient data comprising a set of indices from the 
pulse oximetry readings comprising: delta index, time spent at 
selected oxygen saturation levels and oxygen desaturation events; 

1 0 e) performing random sampling with replacement of patient data 

from (d) to generate a plurality of data sets; 

f) performing multivariate non-parametric analysis in the plurality of 
data sets to generate a plurality of prediction models; 

g) obtaining a provisional predictive AHI value for each patient 
15 corresponding to each prediction model; and 

h) comparing the provisional predictive AHI value to actual AHI 
value for each patient to provide weight to each prediction model, 

wherein the plurality of prediction models and weight of the prediction models 
provides a predictive tool for predicting AHI in an individual. 



20 



2. The method of claim 1 further comprising the step of validating the 
prediction models. 



3. The method of claim 1, wherein the pulse oximetry readings are obtained 
25 overnight. 



4. The method of claim 1 , wherein the selected oxygen saturation is in the 
range of 50 to 100%. 



30 5 The method of claim 4, wherein the selected oxygen saturation levels are 
selected from the group consisting of 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 
92, 93, 94, 95, 96, 97, 98, 99% and combinations thereof. 
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30 



6. The method of claim 5, wherein the selected oxygen saturation levels are 80, 
82, 84, 86, 88 and 90%. 



7. The method of claim 1, wherein the oxygen desaturation levels for 
desaturation events are selected from the group consisting of 2, 3, 4, 5, 6, 7, 8, 9, 10, 
11, 12, 13, 14, 15, 16, 17, 18, 1 9, 20% and combinations thereof. 



8. The method of claim 7, wherein the oxygen desaturation levels are 2, 3 and 
4%. 



9. . The method of claim 1 further comprising: 

a) obtaining additional indices in claim 1(c) selected from the group 
consisting of: electroencephalogram (EEG), electro-oculogram, 
submental and bilateral leg electromyograms, electrocardiogram, 

15 airflow measured by an oral-nasal thermistor, respiratory effort 

measured by thoracoabdominal piezoelectric belts and responses to 
sleep history questionnaires; and 

b) utilizing the additional indices in the step of claim 1(e). 

20 10. A method for predicting AHI in an individual comprising: 

a) obtaining pulse oximetry recordings from the individual; 

b) determining a set of indices from the oximetry readings comprising, 
delta index, time spent at selected oxygen saturation levels and 
oxygen desaturation events; 

25 c) subjecting the set of indices to a plurality of prediction models of 

claim 1(f) to obtain a plurality of provisional predicted AHI values 
corresponding to each model; and 
d) obtaining a final predicted AHI value by averaging the weighted 
provisional predicted AHI values. 



11. The method of claim 1 0, wherein the pulse oximetry readings are obtained 
overnight. 

-19- 



WO 2005/030048 



PCT/US2004/031178 



12. The method of claim 10, wherein the selected oxygen saturation is in the 
range of 50 to 100%. 

13 The method of claim 12, wherein the selected oxygen saturation levels are 
5 selected from the group consisting of 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 
92, 93, 94, 95, 96, 97, 98, 99% and combinations thereof 

14. The method of claim 13, wherein the selected oxygen saturation levels are 
80, 82, 84, 86, 88 and 90%. 

10 

15. The method of claim 10, wherein the oxygen desaturation levels for 
desaturation events are selected from the group consisting of 2, 3, 4, 5, 6, 7, 8, 9, 10, 
11, 12, 13, 14, 15, 16, 17, 18, 19, 20% and combinations thereof. 

15 16. The method of claim 7, wherein the oxygen desaturation levels are 2, 3 and 
4%. 

17. The method of claim 10, further comprising the step of classifying an 
individual as having or not having obstructive sleep apnea based on the final 

20 predicted AHI value. 

1 8. The method of claim 10, wherein step a) further comprises obtaining 
additional indices selected from the group consisting of; electroencephalogram 
(EEG), electro-oculogram, submental and bilateral leg electromyograms, 

25 electrocardiogram, airflow measured by an oral-nasal thermistor, respiratory effort 
measured by thoracoabdominal piezoelectric belts and responses to sleep history 
questionnaires; and step c) further comprises inputting the additional indices into the 
prediction models to obtain the provisional predicted AHI values. 

30 19. A computational system having software thereon to execute the steps 
comprising: 

a) receiving pulse oximetry readings; 
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b) calculating indices comprising: delta index, time spent at selected 
oxygen saturation levels and oxygen desaturation events; 

c) subjecting the indices to a plurality of AHI prediction models to 
generate provisional AHI values corresponding to each prediction 
model wherein each prediction model has specific weight assigned 
thereto; and 

d) calculating a final AHI value based on a weighted average of the 
provisional AHI value. 



10 20. The computational system described in claim 19, wherein: 

a) time spent at select oxygen saturation levels is selected from the 
group consisting of 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 
93, 94, 95, 96, 97, 98, 99% and combinations thereof; and 

b) oxygen desaturation levels for desaturation events are selected from 
15 the group consisting of 2, 3, 4, 5, 6, 7, 8, 9, 10, 1 1, 12, 13, 14, 15, 16, 

17, 18, 19, 20% and combinations thereof. 



2 1 . The computational system of claim 19 further comprising software wherein 
step a) receives additional indices selected from the group consisting of: 

20 electroencephalogram (EEG), electro-oculogram, submental and bilateral leg 

electromyograms, electrocardiogram, airflow measured by an oral-nasal thermistor, 
respiratory effort measured by thoracoabdominal piezoelectric belts and responses to 
sleep history questionnaires; and step c) subjects the additional indices to the 
prediction models. 

25 

22. A computational system having software thereon to execute the steps 
comprising: 

a) receiving indices comprising: delta index, time spent at selected 

t 

oxygen saturation levels and oxygen desaturation events; 
30 b) subjecting the indices to a plurality of AHI prediction models to 

generate provisional AHI values corresponding to each prediction 
model wherein each prediction model has specific weight assigned 
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thereto; and 

c) calculating a final AHI value based on a weighted average of the 
provisional AHI value. 



5 23. The computational system of claim 22, wherein: 

a) time spent at select oxygen saturation levels is selected from the 
group consisting of 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 
93, 94, 95, 96, 97, 98, 99% and combinations thereof; and 

b) oxygen desaturation levels for desaturation events are selected from 
10 the group consisting of 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 

17, 18, 19, 20% and combinations thereof. 



24. The computational system of claim 22 further comprising software wherein 
step a) receives additional indices selected from the group consisting of: 
15 electroencephalogram (EEG), electro-oculogram, submental and bilateral leg 

electromyograms, electrocardiogram, airflow measured by an oral-nasal 
thermistor, respiratory effort measured by thoracoabdominal piezoelectric belts 
and responses to sleep history questionnaires; and step c) subjects the additional 
indices to the prediction models. 

20 
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